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Abstract 

We study the effects of light gluinos with m g ~ 1 GeV on the nucleon parton densities and 
the running of as- It is shown that from the available high-statistics DIS data no lower 
bound on the gluino mass can be derived. Also in the new kinematical region accessible at 
HERA the influence of such light gluinos on structure functions is found to be very small 
and difficult to detect. For use in more direct searches involving final state signatures we 
present a radiative estimate of the gluino distribution in the nucleon. 
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1 Introduction 



The values of the strong coupling constant as(Q) as determined from fixed-target deep 
inelastic lepton-nucleon scattering (DIS) and from hadronic decays of the Z-boson at LEP 
slightly deviate from the running predicted by QCD |l|, |!|. This observation has kindled 
renewed speculations about the possible existence of a light gluino M with a mass of 
at most a few GeV @, ||, 0. The presence of such a light gluino would slow down 
the decrease of ots(Q) between the average scale Q ~ 5 GeV of the DIS measurements 
and the Z-mass and thus improve the overall agreement of experiment with theory ||, [7]. 
The latest a^-determinations have reduced the original discrepancy mainly because of a 
slightly smaller hadronic width of the Z [0. Within errors, the observed running is now 
consistent with standard QCD as well as with the hypothesis of a light majorana gluino. 
Since also other measurement do not convincingly rule out light gluinos ||, this chapter 
is not yet closed. 

If a light gluino exists, it does not only modify the running of as{Q), but also the 
evolution of parton densities and DIS structure functions |], |K]]. It is therefore interesting 
to investigate the question whether or not one can derive lower bounds on the gluino mass 
from present or future structure function measurements alone. Such a bound could not be 
undermined by theoretical uncertainties and would thus put an unavoidable constraint. 
Following ref. |Tl| it is reasonable to concentrate on the mass windows m g ~ 1 GeV 



and 3 GeV ~ m- g ~ 5 GeV. Gluinos in the upper mass range do not significantly change 
the successful description of fixed-target DIS, direct photon and Drell-Yan lepton pair 
production obtained by global fits [12|, |13| in the framework of the QCD improved parton 
model. Furthermore, also extrapolations of DIS structure functions to the new kinematical 
regime accessible at HERA show only little sensitivity to a gluino in this mass range W\ 



However, it is a priori not obvious whether these conclusions also apply to the lower mass 
window m- g < 1 GeV. Such very light gluinos would carry as much as about 10% of the 
nucleon momentum already at scales probed in fixed-target DIS. Thus the task remains 
to clarify whether existing DIS data (including first HERA results at small x fHS| , W 



rule out the existence of light gluinos below 1 GeV. If this is not the case, the question 
arises whether one can expect observable effects later at HERA when more statistics will 
be accumulated in particular at large x and large Q 2 . 

Some of these issues have been investigated in a very recent analysis |l7j following 
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the procedure of the CTEQ collaboration ||13|| . It is concluded that the fixed-target data 



do not discriminate between the presence and absence of light gluinos. On the other 
hand, the scaling violations at x ~ 10~ 2 are found to be distinctly different in these two 
cases, whence small- x data from HERA should soon provide decisive tests. This finding, 



however, contradicts the results of ref. Ul4| for the common mass m~ g = 5 GeV. 



In the present paper, we study the questions raised above within the framework of 
radiatively generating the parton densities from low-scale valence-like input distributions 
|18|| . This approach is especially well suited for the incorporation of a new particle species 



with a mass below the input scale Q Q ~ 2 GeV of conventionally fitted parton distributions 



12] , |13|| , since it avoids the introduction of an additional unknown input function. For 
the gluino, this initial distribution cannot be well determined by a global fit, while our 
procedure leads to a rather definite prediction. This radiative estimate of the gluino 
density in the proton is useful for the study of more direct searches using final state 
signatures || 0, |17| . As far as DIS structure functions are concerned, we find no significant 
effects from gluinos with m g ~ 1 GeV, neither in the kinematical regime explored so far, 
nor in the new region accessible to HERA. Thus our analysis does not confirm the low- a; 



results of ref. 17 mentioned above. 



The paper is organized as follows: In section 2, we recall the formalism for including 
a light gluino in the Q 2 -evolution of hadronic parton densities [0, llOf and derive the 



analytical Mellin-n space solution of the modified evolution equations. In section 3, we 
then explain our phenomenological procedure and present fits to a suitably chosen subset 
of the available high-statistics DIS measurements. The emerging parton distributions 
and their implications are discussed in section 4. Finally, in section 5, we summarize our 
results. 



2 Theoretical Framework 



In this section we formulate the evolution of parton densities in next-to-leading order 
(NLO) in the presence of a light majorana gluino. The corresponding standard QCD 
results can be recovered by simply dropping the gluino-gluon and gluino-quark splitting 
functions together with the gluino effects in the usual splitting functions. All expressions 
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will be given in n-moment space denned by the Mellin transformation 

a n = I dxx n ~ 1 a(x) . 



o 

In this representation, the convolutions in evolution equations and structure functions 
reduce to ordinary products. As a consequence, the equations can be solved analytically. 
The inverse Mellin transformation back to Bjorken-x space has then to be performed 
numerically. 

The strong coupling constant in NLO can be written in the form 

as(Q 2 ) = 1 _ ft lnln(Q 2 /A 2 ) 

4tt ftln(Q 2 /A 2 ) (31 [ln(Q 2 /A 2 )] 2 ' 1 ] 

The (scheme-independent) coefficients of the ftfunction are given by 

A, = 11 ~(2/3)/, ft = 102 - (38/3)/ (3) 

in the standard QCD case, and by 

ft = 9 - (2/3)/ , ft = 54 -(38/3)/ (4) 

in the presence of a majorana gluino. Here / denotes the number of relevant quark 
flavours. The QCD scale parameter A depends on this number as well as on the renor- 
malization scheme. For the present analysis we employ, as usual, the MS scheme. As can 
be seen from eqs. (3) and (4), a light gluino slows down the running of as(Q 2 )- 

Beyond the leading logarithmic approximation parton densities and their evolution 
depend on the factorization scheme, i.e. on the definition of the coefficient functions C, 
entering the moments of the familiar structure functions F\ and F2: 

( ^"(Q 2 ) = f 4 (! + K(^) + ?t(Q 2 )l + ^)^(Q 2 )} , (5) 

where 5i = 2, 5 2 = 1/x. (<&) and G denote the (anti-)quark and gluon distributions, 
respectively, and is the electric charge in units of the electron charge. In contrast to ref. 
[TO] ], we adopt the MS scheme [I9~|. This allows for a direct comparison of the resulting 



parton densities with distributions from conventional QCD analyses. Since the gluino 
coupling to electroweak gauge bosons is the NLO coefficient functions C" are the 



same with and without a gluino. The explicit expressions in MS can be found in ref. |20" |. 



In the MS scheme the massive quark flavours and the gluino can be included in the 
evolution as follows [21]: Below the threshold Q 2 = ml the heavy degree of freedom h is 
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neglected in the evolution of parton densities as well as in the coefficients (3) and (4) of 
the /3-function. At Q 2 > ml it participates in the evolution and the running of cts(Q 2 ) 
like a massless flavour but with the boundary condition 



h n (m 2 h ) = 



(6) 



This guarantees the continuity of the parton densities at threshold. The MS scale pa- 
rameter A^) is changed at Q 2 = m\ by the amount necessary to ensure continuity of 

«s(Q 2 )- 

In the presence of a gluino the evolution equations read 

d .. . > a s (Q 2 ) f p ^)n + as(Q 2 ) p (i)n\ „».(q2) 



cilnQ 2 
d 



r(Q 2 ) 



2tt 
«s(Q 2 ) 



2tt 

p(0)n + ^Sm p{ l)n\ r{Q2) 
2lT J 



, . . , (7) 

d\nQ 2J " 2tt \ 2t ' ' " 

We have simplified the equations by introducing the usual non-singlet combinations of 
quark flavours (i = 1, . . . , f), 



k=l 



together with the vector of singlet moments 



q' 



G n 
\ 9 n 



(9) 



/ being again the number of active quark flavours. The moments of the splitting functions 
are generically denoted by p(°) n and P^ n in first and second order, respectively. Note 
that no additional non-singlet combination involving the gluino distribution g n (Q 2 ) occurs 
since we only consider major ana gluinos here. 

The splitting functions in eq. (7) can be expressed in terms of the standard QCD 
quark and gluon splitting functions by a change of colour factors. In lowest order, P^ n 
is unchanged, while the modified singlet matrix _p(°) n is given by [[J 



p(0)n 



o \ 



p, 



(0)n „(0)n 



V 



Gq 





GG 

6P, 



(0)n 



(10) 



I p (0)n 
l r Gq 

qG 4 J qq / 

Note the absence of quark-gluino splitting in LO and the modification of the gluon-gluon 
splitting function by the gluino contribution to the gluon self-energy graph. The main 
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effect of the gluino on the quarks and gluons is a reduction of the gluon density. The 
second order quantities P^ n and pw™ follow from those without a gluino by changing 
the group factors as described in table 1 of ref. [|IU|Q- The explicit expressions for integer 
values of n can be inferred from the appendix of ref. |2(J using p(°) n = — 'y(°) n /4 and 
p(i)« — —/yCQnys. The analytic continuation to complex n, necessary for the Mellin 



inversion of eq. (1), is explained in ref. [p2 



The solution of the non-singlet part of eq. (7) is straightforward. One finds 



v n {Q 2 ) = (l - a s(<2 2 ) ~ a s(Ql) f p (i )n _ fo_ p (o)n) ] ( as(Q 



7l(3 



2/V 



2 n \ -(2//3 )P< 



(0)n 



111 



In the singlet case, the situation is more complicated due to the non-commutativity of 
the matrices P^ and P^\ An implicit solution is given by 



as(Q 2 )/a s (Ql 
as(Ql 



2tt 



-(2//Sd)p(°)« a s (Q 2 ) 
+ 2tt 

-(2//3b)P<°>" 



«5(q 2 )/« 5 (Qo; 



-(2//3 )p(°>™ 



<xs(Q)/oc s (Qt>) 



U n \q n {Qi) , 



(12) 



where the subleading evolution matrix U is fixed by the commutation relation 



[U n p(°) n l = —jj n -\- pj 1 p n = _p( 1 ) n _ p(°) w 
2 2/5 



)n A 



(13) 



In order to determine {7, we expand the matrix (10) following ref. [E3 



(14) 



i=l 



with A" denoting the eigenvalues of p(°) n . The arbitrarily normalized eigenvectors 

/ 2fP$/(\ i - P®>) ^ 
1 

,(o) 



(15) 



build up the projection matrices 



(16) 



1 There are two misprints in table 1 of ref. pj: the CfNf entries of GL and NSl should both read 
Cx{N F + N x ). 
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where v, 



Substituting U = Ynj=i hVtj and (14) into eq. (13) one readily derives 



A) § " l " % \ - Ai - /3 /2 



(17) 



For simplicity of writing, we have suppressed the superscript n in (15-17). Eqs. (12) and 
(17) then yield the explicit singlet solution 



t \ a s(Qo) 



a S (Ql) ~ a s (Q 2 ) 2 
A) 



a S (Q§) « 5 (Q 2 ) fa s (Q 



2tt 

2 n\ -2(A?-A?)//V 



e?R n e] 



2n 



2tt \a s (Q 



\n-\?-{3 /2 



r(QD ■ 



Finally, the individual quark distributions can be obtained from the singlet and non- 
singlet moments (11) and (18) using the relation 



g?(Q 2 ) + <Z(Q 2 ) = 7 s"(Q 2 ) - t^-iCO 2 ) + E 777— n^-iCO 8 



(19) 



3 Phenomenological Analysis 

We begin with the question whether or not a very light gluino is consistent with the high- 
statistics fixed-target DIS measurements. To find the answer it is not necessary to carry 
out a complete global fit to all data available. Instead we follow ref. and proceed in 
two steps. First, we study the valence quark region and the strong coupling constant. 
Using the result as an input we then investigate the more interesting regime at smaller x 
where the gluon and sea distributions play an important role. Only this second part of 
the analysis requires new fits to data. 

For the valence quark description of F<i and xF^ in standard QCD we adopt the NLO 
proton distributions and the QCD scale parameter from the global fit KMRS(B_) p4[fl. 
At the reference scale Ql = A GeV 2 , they are given by 

x(u v + d v )(x,Ql) = 0.385 s a27 (l + 9.9^+ 17. 7x)(l-x) 3 - 93 

xd v {x,Q 2 ) = 1.50x a61 (l + 1.08v^)(l -x) 4 - 68 ( 20) 



2 In refs. p4| and |2?| a slightly different NLO relation between as and A has been used. We have 
transformed their results to the convention (2) employed here, resulting in a 10 MeV increase of A< 4 >. 
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and A^ 4 ) = 200 MeV. More recent fits exist, the differences in the valence quarks being, 
however, marginal and completely immaterial for the present purpose. Also the value of 
A^ is consistent with more recent determinations, for example from newer CCFR data 
giving A^ 4 ) = (220 ±50) MeV (2^| 2 . With a gluino in the mass range m~ g < Qq included it 
is still possible to reproduce xF 3 (x,Q 2 ) and the valence contribution to F 2 (x, Q 2 ) in the 
measured range, 5 GeV 2 £ Q 2 Z 300 GeV 2 and 0.01 £ x Z 0.7, within less than 0.5 % by 
a slight modification of the valence quark input and a suitable choice of A~ 4 \ Hence no 
new valence fit to data is necessary. Optimization of this agreement leads to 

x{u v + d v ){x,Ql)- g = 0.388a; a27 (l + 9.9v / x + 17.6a;) (1-x) 3 - 975 

xd v (x,Ql)~ g = 1.51x a61 (l + 1.09>/x)(l -a;) 4 ' 74 (21) 

together with Af = 30 MeV. 

The strong coupling constants ats{Q) i n the two scenarios (20) and (21) coincide at 
Q 2 = 12 GeV 2 . The continuity of as at the heavy quark thresholds Q 2 = m 2 (h = c, b) 
then implies for the NLO scale parameters A^ in the MS scheme: 

A (3,4,5) = 24g; 200; 131 Me y 

A j3,4,5) = 5g 1Q g ^ ^ 

where m c = 1.5 GeV, m& = 4.5 GeV and m g < m c has been used. Obviously, these values 
of A^ are independent of the gluino mass, but apply only above the gluino threshold 
Q = m,g. The gluino mass dependent values of A^ relevant below the gluino threshold 
will be given later. Using eq. (22) we find as(M z ) = 0.128 (0.109) with (without) the 
gluino. The result in brackets agrees with the recent CCFR value, as(Mz) = 0.111±0.005 



25(1 obtained from a standard QCD analysis. A similar uncertainty in ots(Mz) is expected 
for the gluino case. The existence of a light gluino would, however, not only modify the 
DIS results on as(Q) and the extrapolation to Q = Mz but also change the direct LEP 
determinations of as(Mz) due to contributions to virtual and real radiative corrections 
0. With (without) a light gluino one obtains §, a s {M z ) = 0.132 (0.120) ± 0.006 from 
hadronic event shapes and as(Mz) = 0.124 (0.122) ± 0.009 from the hadronic width r\ 
of the Z°. Note that the original results of ref. have been updated with respect to the 
latest value of IV Cleary, within present uncertainties, both scenarios give a consistent 
picture. 

We now turn to the more interesting part of the analysis dealing with the gluon and 



sea densities. For this, we adopt the radiative approach developed in refs. p2|, ITR |26|. [27[ . 
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This approach provides a unified description of hadronic and photonic parton distributions 
which agrees with experiment, and which predicted the steep rise of the structure function 
F 2 at small x discovered recently at HERA [T5|, [TJJ. The charactistic features of this 
framework are a very low scale (///A®) 2 ~ 5 -=-6 at which the Altarelli-Parisi evolution is 
started and valence-like input densities. For light gluinos with m~ g > fi, this procedure has 
the advantage that the gluino can approximately be treated like a massive quark with the 
boundary condition (6). It is generated similarly to the strange sea discussed below. In 
this way one can avoid the introduction of a gluino input distribution except for almost 
massless gluinos with rrig < fi. For the gluon and sea quark inputs we use |H 



xG(x,fi 2 ) = Ax a {l-xf 

xu{x,n 2 ) = xd{x,fi 2 ) = A'x a '{l -xf (23) 
xs(x, fj, ) = xs(x, /i 2 ) = . 

The vanishing of the strange sea at /i 2 leads to SU(3) breaking in the sea in accordance 
with experimental results from opposite-sign dimuon events in neutrino-nucleon scattering 
|28fl . The value of /i is fixed by the momentum fraction <x> v (/i 2 ) carried by the valence 



quarks. For the standard QCD case we take fx 2 = 0.3 GeV 2 corresponding to <x> v ~ 0.6 
fl8fl . Keeping this fraction also in the presence of a gluino, the scale /i decreases with 
the gluino mass due to the slower evolution at low Q 2 . This correlation is shown in table 
1 for typical gluino masses. In order to illustrate the sensitivity of the outcome of the 
subsequent fits to the exact value of \i we also give results for /i 2 = 0.36 GeV 2 , that is, 
for a 10% increase of /i. 

The input parameters appearing in (23) are fitted to the available high-statistics data 
on Ff and from SLAC |2jJ, BCDMS [00 and NMC |2[ in the region of good 



sensitivity, i.e. at x < 0.3. Statistical and systematical errors are added in quadrature. 

As in ref. [|T2] the BCDMS data is normalized down by 2% relative to the SLAC results. 

In order to suppress possible higher twist effects we restrict our fits to measurements at 

Q 2 > 5 GeV 2 . An additional cut on the invariant mass W of the hadronic final state 

is not neccessary in view of the low a;-values involved. In total, we are left with 147 

data points. Additional high precision data on F^ N , F% and the ratio F^/F^ mainly 

influence the flavour decomposition, and are unimportant for the gluino issue. Moreover, 

by using the valence distributions from ref. p4j the main features of these measurements 

are incorporated. Likewise, data on Drell-Yan lepton pair production mainly constrain 

3 We take the BCDMS data as used in |3l]], i.e. with a shift of the central values due to the main 
systematic error. However, in the x-region considered here the effect is not essential. 
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the shape of the sea quark densities at large x where they are small and not relevant for 
the present study. On the other hand, data on direct photon production in proton-proton 
collisions may be relevant, since this process directly probes the gluon distribution 
which is most affected by light gluinos. However, unlike in the case of DIS and the 
Drell-Yan process, the hard photon production cross sections in NLO do get modified in 
the presence of gluinos. These modifications have not yet been calculated. Hence, we 
disregard the direct photon data in the fits for theoretical consistency but take them into 
account approximately as discussed below. 



m~ g 


A< 3 ) /MeV 


^/GeV 2 


<X> V (/J 2 ) 


a' 




A 1 


x 2 


0.4 


118 


0.09 


0.617 


0.490 


8.03 


0.138 


140 


0.7 


148 


0.125 


0.618 


0.519 


8.02 


0.134 


139 


1.0 


168 


0.155 


0.615 


0.514 


8.03 


0.131 


137 


1.3 


184 


0.18 


0.615 


0.525 


8.06 


0.130 


137 


no g 


248 


0.30 


0.614 


0.529 


7.93 


0.127 


134 


0.4 


118 


0.12 


0.572 


0.362 


7.66 


0.134 


134 


0.7 


148 


0.155 


0.578 


0.371 


7.63 


0.130 


134 


1.0 


168 


0.19 


0.576 


0.364 


7.59 


0.129 


133 


1.3 


184 


0.22 


0.575 


0.364 


7.56 


0.128 


132 


no g 


248 


0.36 


0.574 


0.367 


7.44 


0.125 


130 



Table 1: Fits of DIS structure functions to experimental data on F% at x < 0.3 and 
Q 2 > 5 GeV 2 (147 data points) in standard QCD (no g) and including a light gluino. The 
input distributions are given in eq. (23) and specified further in the text. 



Since the gluon normalization A is fixed by the momentum sum rule, we have five free 
parameters describing the input distributions (23). These are fitted to the data selected 
above in the standard QCD case and for several gluino masses m g between 0.4 GeV and 
1.3 GeV. The results of the fits are given in table 1. The outcome is almost insensitive to 
variations of the two shape variables of the gluon density in the range of values favoured 



by direct photon measurements |33|, 34] . For this reason, we simply take 



a = 2.0 (1.5) , $ = 4.0 (24) 

in the fits summarized in the upper (lower) part of the table. Also the sea quark shape 
parameters a' and /3' are only slightly affected by the inclusion of a light gluino. However, 
the overall normalization A' shows a small but significant rise with decreasing gluino 
mass. This is exactly what one expects physically since the gluon density, and hence the 
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evolution of quark-antiquark pairs, is somewhat diminished due to gluino radiation. The 
quality of the above fits is as good as the one of the recent MRS global analysis [fl2 1 , if 
one compares the description of the data used in both fits. As can be seen from the \ 2 
values in the last row of table 1 there is no evidence from the fits in favour or against 



gluinos with masses below 1 GeV. This conclusion agrees with ref. [17 



4 Resulting Parton Distributions 

Using the results from the fits outlined in the previous section, we now discuss the parton 
densities of the nucleon and the structure function F 2 P . This will also enable us to answer 
the second question raised in the introduction concerning the potential of future DIS 
measurements at HERA in testing the existence of a very light gluino. As reference cases 
for comparison we focus on the standard QCD fit with fi 2 = 0.30 GeV 2 and the gluino 
scenario with m g = 0.7 GeV and /i 2 = 0.125 GeV 2 . 

It is instructive to consider briefly the fractions <x>j= ^dxxf{x, Q 2 ) of the nucleon 
momentum carried by quarks, gluons and gluinos. Their Q 2 -evolution is shown in fig. 1. 
As the main effect, the average gluon momentum is lowered by gluino radiation. This 
reduction is fed back to the quarks at large Q 2 . Below threshold, at Q 2 ~ m~, the change 
of the quark and gluon momentum fractions reflects the slower running of ots{Q 2 ) in the 
presence of a gluino. Asymptotically, the gluino momentum reaches about 10% of the 
proton momentum. Moreover, one can see that already at Q 2 ~ 10 GeV 2 the gluino 
component gives a non-negligible contribution to the momentum sum rule. 

Figs. 2 and 3 exhibit the sea quark distributions xu = xd and the gluon density xG 
at Q 2 = 10 GeV 2 . While the inclusion of the gluino has virtually no effect on the sea 
quarks, the gluon distribution is reduced by about 15% over the range x ~ 0.05. This 
shrinkage is, however, not expected to lead to phenomenological problems, most notably, 
in direct photon production, given the present theoretical and experimental uncertainties 
and taking into account the partial compensation by the larger value of ots{Q 2 ) at high 

Q 2 - 

The predictions in the small-x region and the Q 2 -evolution of xu and xG are depicted 
in figs. 4 and 5. Also here, the gluino effect on the sea quark density can be neglegted 
for all practical purposes. Furthermore, the depletion of the gluon distribution shown in 
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fig. 3 for large x disappears as x decreases. All effects depend only weakly on Q 2 . In 
particular, in the region of x and Q 2 currently under investigation at HERA, it appears 
impossible to uncover a light gluino in the structure functions. This conclusion is further 
corroborated in fig. 6, where is displayed in this kinematical region together with 



first HERA data [ 15| , [16| . The tiny down-shifts of F$ are mainly due to a smaller charm 
contribution^ which in turn is a consequence of the lower gluon density. The standard 
QCD results with fx 2 = 0.3 GeV 2 are practically identical to the q pre-HERA predictions 



of refs. [pl|. The uncertainty in these predictions is illustrated by the expectations for 



H 2 = 0.36 GeV 2 . 



The above findings differ from the results obtained in ref. JIT]. While, according to 
this paper, at Q 2 = 10 GeV 2 the quark and gluon distributions are only affected at the 
10% level, at Q 2 = 100 GeV 2 and 10~ 4 ~ x ~ 10 -3 the presence of a light gluino gives rise 
to an increase of xu and xG by about 50%. Moreover, this strong evolution is roughly 



independent of the gluino mass for m~ g ~ 5 GeV. This disagrees also with ref. W% -, where 
a gluino of 5 GeV is shown to change F 2 P by at most a few per cent. 

Turning to the Q 2 -evolution of structure functions at large x, in fig. 7 we illustrate the 
changes in F 2 P in the region Q 2 /x < s ~ 10 5 GeV 2 which will be probed at HERA. The 
deviations from the standard QCD evolution are somewhat bigger than the effects found 



in ref. for rrig = 5 GeV. This is expected because of the lower gluino mass considered 



here. However, also for such light gluinos the deviations do not exceed a few per cent and 



will therefore be very hard to detect ||36|| , in contrast to the expectation in ref. J7J based 
on the enhancemant of as at high scales. This conclusion suggests to return once more to 
the gluino content in the nucleon considered at the beginning of this section which may 
play a role in more direct searches at ep- and hadron colliders. 

Fig. 8 shows distributions in x and the evolution of the gluino component together 

with the strange sea also generated radiatively from a vanishing input. The threshold 

dependence is illustrated for the gluino by considering three different masses m g = 0.4, 

0.7 and 1.0 GeV. At large x, g(x,Q 2 ) is about three times bigger than s(x, Q 2 ) and, 

incidentally, rather similar to u(x, Q 2 ). The ratio g/s rises with decreasing x, reaching 

a value around five at x ~ 10 -4 in accordance with the expectation from the differences 

in the colour factors of the relevant splitting functions in eq. (10). The relative abun- 

4 Above threshold the charm quark is treated as a massless flavour. Using instead the complete photon- 
gluon fusion cross section for the charm contribution results in a decrease of F% by 5 -j- 10 % at small x. 



For a detailed discussion see |35 
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dance of gluinos depends only weakly on Q 2 except close to the threshold. Comparing 
these predictions with the gluon and gluino distributions given in ref. [|17]] we find rough 
agreement at Q 2 = 10 GeV 2 . However, at Q 2 = 100 GeV 2 and small x a gluino with 
rrig ~ is predicted to be about three times as abundant as what we find. Moreover, it 
dominates the corresponding up-quark density by almost a factor of nine in contrast to 
the expectation from eq. (10). We have checked that our evolution program reproduces 
the results of ref. [Q if the same input is used. Therefore, we believe that fig. 8 represents 
a reasonable estimate of the gluino content of the proton, if a light gluino exists. 



5 Summary 

We have investigated the effects of a very light gluino with m g ~ 1 GeV on the quark and 
gluon distributions of the nucleon and the DIS structure functions. The whole range in 
Bjorken-x and Q 2 which is being probed in fixed-target experiments and at HERA has 
been taken into consideration. We have adopted the framework described in refs. p2| , |T8] 



in which the parton densities are generated radiatively from valence-like distributions 
at a very low input scale. This procedure leads to rather definite predictions at small 



x which are in agreement with first HERA results ^5j, [16| . Moreover, in the present 
application it has the advantage that the light gluino can be included in the Q 2 -evolution 
similarly as a massive quark, i.e. without introducing a perturbatively uncalculable and 
phenomenologically rather unconstrained gluino input function. 

( f) 

Using as determined from a valence quark analysis of F 2 and xF 3 in the presence 
of a light gluino, the valence-like input distributions have been fitted to the relevant high- 
statistics data on F£- Globally, we find that the gluino carries about 5 -f- 10 % of the 
nucleon momentum. On the other hand, the effects on the quark and gluon distributions 
at x Z 1Q- 2 are completely negligible, even after evolution over a few orders of magnitude 
in Q 2 . The fits to fixed-target data including a light gluino are almost as good as the 
standard QCD fits. Moreover, we have shown that the structure function F 2 P is practically 
insensitive to light gluinos also in the new kinematical region covered by HERA. We 
thus conclude that measurements of DIS structure functions are not able to discriminate 
between the presence and absence of light gluinos. 

Then, in order to test the light-gluino hypothesis one has to study final state signatures 
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such as jet rates and and displaced vertices. In pp, pp and ep collisions, this involves the 
gluino density in the proton. We have presented a radiative estimate for this component 
assuming that the gluino distribution vanishes at the threshold Q 2 = m~. Far from 
threshold the gluino density is predicted to be about three to five times higher than the 
strange quark density. 
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Figure 1: The NLO Q 2 -evolution of the momentum fractions <x>f carried by quarks 
and gluons in standard QCD and including a light gluino with m~ g = 0.7 GeV. The gluino 
momentum fraction is given for varying gluino masses. 
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Figure 2: NLO sea quark distribution with and without a light gluino in the x-region 
probed by fixed-target DIS. 
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Figure 6: NLO predictions for the proton structure function F 2 P at small x in comparison 
to first HERA data from the ZEUS |15] and HI |16|] collaborations. The solid curves 
represent the standard QCD results for the input scales /i 2 = 0.3 GeV 2 (upper) and 
yU 2 = 0.36 GeV 2 (lower). The dashed curves show the corresponding predictions with a 
light gluino included. 
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Figure 7: The ratio of the structure functions F 2 P with and without a light gluino as 
extrapolated from our fits into the large-rr, large-Q 2 region accessible at HERA. 
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Figure 8: The radiatively generated NLO strange quark and gluino distributions of the 
proton. The gluino effects on xs are shown for a fixed gluino mass, while xg is given for 
varying gluino masses. 
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Figure 8: The radiatively generated NLO strange quark and gluino distributions of the 
proton. The gluino effects on xs are shown for a fixed gluino mass, while xg is given for 
varying gluino masses. 
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